Influence of frozen capillary waves on contact mechanics 
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Free surfaces of liquids exhibit thermally excited (capillary) surface waves. We show that the 
surface roughness which results from capillary waves when a glassy material is cooled below the 
glass transition temperature can have a large influence on the contact mechanics between the solids. 
The theory suggest a new explanation for puzzling experimental results [L. Bureau, T. Baumberger 
and C. Caroli, arXiv:cond-mat/0510232 about the dependence of the frictional shear stress on the 
load for contact between a glassy polymer lens and flat substrates. It also lend support for a recently 
developed contact mechanics theory. 



Many technological applications require surfaces of 
solids to be extremely smooth. For example, window 
glass has to be so smooth that no (or negligible) diffu- 
sive scattering of the light occur (a glass surface with 
strong roughness on the length scale of the light wave- 
length will appear white and non-transparent because of 
diffusive light scattering). For glassy materials, e.g., sili- 
cate glasses, or glassy polymers, e.g., Plexiglas, extremely 
flat surfaces can be prepared by cooling the liquid from 
well above the glass transition temperature T g , since in 
the liquid state the surface tension tends to eliminate (or 
reduce) short-wavelength roughness. Thus, float glass 
(e.g., window glass) is produced by letting melted glass 
flow into a bath of molten tin in a continuous ribbon. 
The glass and the tin do not mix and the contact sur- 
face between these two materials is (nearly) perfectly flat. 
Cooling of the melt below the glass transition tempera- 
ture produces a solid glass with extremely flat surfaces. 
Similarly, spherical particles with very smooth surfaces 
can be prepared by cooling liquid drops of glassy mate- 
rials below the glass transition temperature. Sometimes 
glassy objects are "fire polished" by exposing them to a 
flame or an intense laser beamQ which melt a thin sur- 
face layer of the material, which will flow and form a very 
smooth surface in order to minimize the surface free en- 
ergy. In this way, small-wavelength roughness is reduced 
while the overall (macroscopic) shape of the solid object 
is unchanged. 

However, surfaces of glassy materials prepared by cool- 
ing the liquid from a temperature above T g cannot be 
perfectly (molecularly) smooth, but exhibit a fundamen- 
tal minimum surface roughness, with a maximum height 
fluctuation amplitude of typically 10 nm, which cannot 
be eliminated by any changes in the cooling procedure. 
The reason is that at the surface of a liquid, fluctua- 
tions of vertical displacement are caused by thermally 
excited capillary waves (ripplons)jll3. At the surface of 
very viscous supercooled liquids near the glass transition 
temperature these fluctuations become very slow and are 
finally frozen in at the glass transition^- 

The roughness derived from the frozen capillary waves 
is unimportant in many practical applications, e.g., in 



most optical applications as is vividly evident for glass 
windows. However, in other applications they may be of 
profound importance. In this letter I show that they are 
of crucial importance in contact mechanics. For elastic 
hard solids such as silica glass, already a roughness of 
order 1 nm is enough to eliminate the adhesion, and re- 
duce the contact area between two such surfaces to just a 
small fraction of the nominal contact area even for high 
squeezing pressures 0. In this letter I will show that this 
is the case even for elastically much more compliant solids 
such as Plexiglas (PMMA). The results presented below 
suggest a new explanation for the puzzling experimen- 
tal result of Bureau et al0 for PMMA, and in addition 
lends support for a recently developed contact mechanics 
theory^]. Bureau et al observed that the nominal shear 
stress <7 when a PMMA lens was slid on a flat substrate, 
depends strongly on the nominal pressure (or perpendicu- 
lar stress) Pnom- They assumed that perfect (molecular) 
contact occurred in the nominal contact area Aq, and 
interpreted the experimental results as a strong depen- 
dence of the frictional shear stress on the local pressure. 
However, the analysis below shows that the experimental 
results can be understood assuming that the area of real 
contact A is smaller than the nominal contact area Ao, 
and assuming a constant (pressure independent) shear 
stress <jf in the area of real contact in such away that 
the friction force Ff = at A = uAq and the normal load 
Fn = pA — Pnom^Oj where p is the (average) pressure in 
the area of real contact. 

The most important property of a rough surface is the 
surface roughness power spectrum which is the Fourier 
transform of the height-height correlation function^]: 

c(q) = l^f\ d2x W x )M°)> e " q ' x 

Here z = h(x) is the height of the surface at the point 
x = (x, y) above a flat reference plane chosen so that 
(h(x)) = 0. The angular bracket (...) stands for ensemble 
averaging. The surface roughness power spectra due to 
capillary waves is of the form0, 0, Q 

C< " q) = (2^)2 pg + lq 2 + Kq 4 ^ 



where 7 is the surface tension, k the bending stiffness, 
and p the mass density of the glassy melt. The smallest 
relevant wavevector q = qq « 2n/L, where L is the diam- 
eter of the nominal contact area Aq. In the applications 
below L « 0.1 mm and (for PMMA) 7 w 0.04 J/m 2 , 
giving 7g 2 > 10 6 J/m 4 >> pg. Thus, the gravity term 
in (1) can be neglected. The mean of the square of the 
surface height fluctuation is given by 
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where q c = (j/n) 1 / 2 is a cross-over wavevector. For the 
polymer surface which interests us below (PMMA), Eq. 
(2) gives the rms roughness ss 0.5 nm which correspond 
to a maximum roughness amplitude of about 3 nm. I 
will now show that even this small roughness has a large 
influence on the contact mechanics. 




can also be applied when the adhesional interaction is 
included. 

Fig-Dshows the contact between two solids at increas- 
ing magnification £. At low magnification (£ = 1) it 
looks as if complete contact occurs between the solids at 
many macro asperity contact regions, but when the mag- 
nification is increased smaller length scale roughness is 
detected, and it is observed that only partial contact oc- 
curs at the asperities. In fact, if there would be no short 
distance cut-off, the true contact area would vanish. In 
reality, however, a short distance cut-off will always exist 
since the shortest possible length is an atomic distance. 
In many cases the local pressure at asperity contact re- 
gions at high magnification will become so high that the 
material yields plastically before reaching the atomic di- 
mension. In these cases the size of the real contact area 
will be determined mainly by the yield stress of the solid. 

The stress distribution P(a, £) at the interface when 
the interface is studied at the magnification £ = L/X, 
where L is the diameter of the nominal contact area be- 
tween the solids and A the shortest surface roughness 
wavelength which can be detected at the resolution £, 
satisfies the diffusion-like equation 



where 



(4) 
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where E is Young modulus and v the Poisson ratio, and 
where qo — 2ir/L and q — £go- The stress distribution 
function P(cr, £) must vanish for a < — cr a , where <r a (£) is 
the highest tensile stress possible at the interface when 
the system is studied at the magnification £. Substituting 
() in () gives 



FIG. 1: An elastic (e.g., rubber) block (dotted area) in adhe- 
sive contact with a hard rough substrate (dashed area). The 
substrate has roughness on many different length scales and 
the rubber makes partial contact with the substrate on all 
length scales. When a contact area is studied at low magnifi- 
cation (£ = I) it appears as if complete contact occurs in the 
macro asperity contact regions, but when the magnification 
is increased it is observed that in reality only partial contact 



Recently, a contact mechanics theory has been devel- 
oped which is valid not only when the area of real contact 
is small compared to the nominal contact area, but which 
is particularly accurate when the squeezing force is so 
high that nearly complete contact occurs within the nom- 
inal contact areaM lift 1 1 1 1 | . All other contact mechanics 
theories [H El Ej were developed for the case where 
the area of real contact is much smaller than the nomi- 
nal contact area. The theory developed in Ref. 0, 0] 
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If one assume, as is expected from the theory of cracks, 
that the (tensile) stress diverge as ~ r -1 / 2 with the dis- 
tance r from the boundary line between a contact island 
and the surrounding non-contact region, then one can 
easily prove the boundary condition P{— c a (C),C) = 0- 
When the adhesion is neglected cr a = and P(0, Q = 0. 
The (normalized) area of (apparent) contact when the 
system is studied at the resolution £ is given by 
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where Aq is the nominal contact area. The area of real 
contact is obtained from (7) at the highest (atomic) res- 
olution, corresponding to the magnification £1 = L/a, 
where a is an atomic distance. 



The detachment stress er a can be determined using the 
theory of cracks which gives 

, n „ ( je S (C)Eq \ 1/2 

CTa(c) *(,27T^)J () 

Here 7 o ff(0 is the (effective) interfacial binding energy 
when the system is studied at the magnification (. At 
the highest magnification (corresponding to atomic reso- 
lution) 7 c fj = A7, where A7 = 71 +72—712 is the change 
in the energy (per unit area) when two flat surfaces of the 
solids are brought into adhesive contact. The effective in- 
terfacial energy at the magnification £ is determined by 
the equation 

A(C)7eff(C) = ^(Ci)A7-f4i(C) (9) 

where C/ e i(C) is the elastic energy stored at the interface 
due to the elastic deformation of the solid on length scales 
shorter than A = L/Q, necessary in order to bring the 
solids into adhesive contact. An explicit expression for 
7eff is given in Ref. [Tof . 

In the experiment by Bureau et al.Q, the PMMA lens 
were prepared by cooling a liquid drop of PMMA from 
250 °C to room temperature. In the liquid state the sur- 
face fluctuations of vertical displacement are caused by 
thermally excited capillary waves (ripplons). When the 
liquid is near the glass transition temperature T g (about 
100 °C for PMMA), these fluctuations become very slow 
and are finally frozen in at the glass transition. Thus, 
the temperature T in (6) is not the temperature where 
the experiment was performed (room temperature), but 
rather the glass transition temperature T g . 

The (upper) solid line in Fig. [21 shows the calculated 
relative contact area as a function of the (nominal) pres- 
sure Pnom = Fn/A (where i*N is the squeezing force). In 
the calculation we have used the equations above with 
the measured elastic (Young) modulus E = 2.9 GPa 
and surface tension 7 = 0.04 J/m 2 . We have also 
used the glass transition temperature T g w 370 K, the 
short-distance cut-off wavevector qi = 7 x 10 9 m -1 and 
q c = 1.7 x I0 9 m -1 . The cross-over wavevector q c (or the 
bending stiffness n) has not been measured for PMMA, 
but has been measured for other systems. Thus, for alka- 
nes at T « 100 °C for C20 and C36, q c « 4.4 x 10 9 and 
« 2.7 x 10 9 m _1 , respectively 1 15| . There are also some 
studies of polymers [lflUJi (using the atomic force micro- 
scope) giving ijc k 1 x 10 8 m _1 . The interfacial binding 
energy A7 can be estimated using0 A7 w 2(71 72) 1 / 2 ~ 
0.07 J/m 2 , where 71 = 7 w 0.04 J/m 2 is the surface en- 
ergy of PMMA and 72 « 0.02 J/m 2 the surface energy of 
the (passivated) substrate. In the calculations we used 
A7 = 0.06 J/m 2 . 

If one assumes that, in the relevant pressure range (see 
below), the frictional shear stress Cf between PMMA and 
the substrate is essentially independent of the normal 
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nominal pressure (MPa) 

FIG. 2: The calculated relative contact area as a function of 
the (nominal) pressure p n om including adhesion (upper solid 
line) and without adhesion. Circles: the normalized, (nomi- 
nal) shear stress a /at for PMMA sliding on TMS as a func- 
tion of the (nominal) pressure. The nominal shear stress a 
has been divided by 01 = 50 MPa, which is the (measured) 
shear stress in the area of real contact (see text). Squares: the 
normalized, (nominal) shear stress a /at (with at = 5 MPa) 
for PMMA sliding on OTS as a function of the (nominal) 
pressure. 

stress p in the asperity contact regions, then the calcu- 
lated curve in Fig. 121 is also the ratio between the nom- 
inal (or apparent) shear stress a and the true stress at; 
A/Aq = cr/cri (note: the friction force Ff = <jfA = aAg). 

The circles and squares in Fig. [3 show the measured 
data of Bureau et al.[j| for er/o-f. They performed ex- 
periments where a PMMA lens, prepared by cooling a 
liquid drop of PMMA, was slid on silicon wafers (which 
are nearly atomically smooth) covered by a grafted silane 
layer. Two different types of alkylsilanes were employed 
for surface modification, namely a trimethylsilane (TMS) 
and octadecyltrichlorosilane (OTS). The circles in Fig. |3 
are the shear stress, divided by <jf = 50 MPa, for PMMA 
sliding on TMS as a function of the (nominal) pressure 
Pnom- The shear stress at = 50 MPa was deduced from 
multi-contact experiments^^ (where the local pressure 
in the asperity regions is so high as to give rise to lo- 
cal plastic deformation) using Of w crH/p uom , where 
H w 300 MPa is the hardness (yield stress) of PMMA as 
determined from indentation experiments. This equation 
follows from oAq = a^A and AH = AQp nom . The squares 
in Fig. Ushow a/a f for PMMA sliding on OTS as a func- 
tion of the (nominal) pressure p n0 m- For this system no 
multicontact experiments were performed, and we have 
divided the apparent shear stress by 01 — 5 MPa, chosen 
so as to obtain the best agreement between the theory 
and experiment. Hence, in this case at represent a theo- 
retically predicted shear stress for the multicontact case; 
we suggest that multicontact experiments for the system 
PMMA/OTS are performed in order to test this predic- 



tion. 

The calculation presented above shows that the lo- 
cal pressure in the asperity contact regions is of or- 
der ~ 2 GPa, which is much higher than the macro- 
scopic yield stress (about 300 MPa) of PMMA. However, 
nanoscale indentation experiments |2f| have shown that 
on the length scale (~ 10 nm) and indentation depth 
scale (~ 1 nm) which interest us here, the indentation 
hardness of PMMA is of order ~ 10 GPa, and no plastic 
yielding is expected to occur in our application. 

The study above has assumed that the frictional shear 
stress at the PMMA-substrate interface is independent 
of the local pressure p, which varies in the range p 
50 — 100 MPa. This result is expected because the 
shear stress will in general exhibit a negligible pressure 
dependence as long as p is much smaller than the ad- 
hesional pressure p a( j. The adhesional pressure is de- 
fined as follows: In order for local slip to occur at the 
interface, the interfacial molecules must pass over an 
energetic barrier of magnitude Se. At the same time 
the spacing between the surfaces must locally increase 
with a small amount (some fraction of an atomic dis- 
tance) Sh. This correspond to a pressure work pa Sh, 
where a 2 is the contact area between the molecule, or 
molecular segment, and the substrate. Thus the total 
barrier is of order Se + pa 2 Sh = a 2 5h{p &< \ + p) where 
Pad = Se/(a 2 Sh). For weakly interacting systems one 
typically have[2lj Se/a 2 « 1 meV/A 2 and Sh « 0.01 A 
giving p a d « 10 GPa, which is at least one order of magni- 
tude larger than the pressures in the present experiment. 
We note that the pressure p a d is similar to the pressure p 
which must be applied before the viscosity of liquid hy- 
drocarbon oils start to depend on the applied pressure. 
The reason for the pressure dependence of the viscosity 
of bulk liquids is similar to the problem above, involving 
the formation of some local "free volume" (and the as- 
sociated work against the applied pressure) in order for 
the molecules to be able to rearrange during shear. 

To summarize, I have shown that the surface roughness 
which result from capillary waves when a glassy material 
is cooled below the glass transition temperature can have 
a large influence on the contact mechanics between the 
solids, even for relatively (elastically) soft solids such as 
PMMA. This fact suggest a new explanation for puz- 
zling experimental results about the dependence of the 
frictional shear stress on the load for contact between a 
PMMA lens and flat passivated silica substrates. It also 
lend support for a recently developed contact mechanics 
theory 
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If one assumes that the energy necessary to move over 
the barrier 8e is lost during each local slip event, then the 
shear stress o~i ~ (5e/a 2 )/b where b is the slip distance 
which will be of order a lattice constant. Thus, from the 
measured at (about 50 MPa for TMS) and using 6 sa 
3 A gives Se/a 2 w a { b » 1 meV/A 2 for TMS. For OTS 



8e/a 2 and the adhesion pressure would be about 1 order 
of magnitude smaller than for TMS. This analysis does 
not take into account that the (measured) shear stress 
erf involves thermally- assisted stress-induced dissipative 



events (i.e., at would be larger at lower temperatures), 
so the estimation of Se/a is a lower bound. 



